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Abstract 

> 

\& • An exclusion particle model is considered as a highly simplified 

model of a limit order market. Its price behavior reproduces the well 
\^ | known crossover from over-diffusion (Hurst exponent H > 1/2) to dif- 

<^ ' fusion {H = 1/2) when the time horizon is increased, provided that 

orders are allowed to be canceled. For early times a mapping to the 
totally asymmetric exclusion process yields the exact result H = 2/3 
which is in good agreement with empirical data. The underlying uni- 
versality class of the exclusion process suggests some robustness of the 
exponent with respect to changes in the trading rules. In the crossover 
£5 \ regime the Hurst plot has a scaling property where the bulk deposi- 

tion/cancellation rate is the critical parameter. Analytical results are 
fully supported by numerical simulations. 

PACS numbers: 05.40; 89.90 +n 
Keywords: Limit order market, Hurst exponent, KPZ equation, asymmetric 
exclusion process. 

1 Introduction 

Financial markets have in recent years been at the focus of physicists at- 
tempts to apply existing knowledge from statistical mechanics to economic 
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problems [I], 0, [3], ||. These markets, though largely varying in details of 
trading rules and traded goods, are characterized by some generic features 
of their financial time series, called 'stylized facts' [T|, ^|, ffl. Agent based 
models of financial markets are successful to reproduce some stylized facts 
B H B H) H) H3 ' sucn as volatility clustering, fat-tailed probability distribu- 
tion of price increments and over-diffusive price behavior at short time scales 
and diffusive behavior at later times. But all of them need an explicit price 
formation rule that links excess demand to price changes [11], ||, [L^, |13"1 , that 
can be itself problematic. An other approach consists in modeling the price 
formation, for instance in limit order markets (TTJ], |T5, [IT], [TB[ . So far, all 
these models of limit order markets have under-diffusive prices at short times, 
with a crossover to diffusive prices at longer times for some of them. Under- 
diffusive behavior at short times is realistic in limit order markets, but all 
these models lack the over-diffusive price behavior observed in real markets 
at medium time scales. Here, we introduce a crude non-equilibrium model 
with overdiffusive price that is able to reproduce the crossover from a Hurst 
exponent H = 2/3 to H = 1/2 at larger times, when correlations in the price 
dynamics are washed out by cancellations of existing orders and independent 
placements of new orders. In the early time regime our model belongs to the 
lrf-KPZ universality class [pi|, hence, its mechanism for over-diffusive price 
spreading is robust and analytically tractable. 

In section II we define our model in terms of the limit order market dy- 
namics. Our simulation results are presented in section III. In section IV the 
equivalence of the early-time regime of our model to the totally asymmetric 
exclusion process (TASEP) p0| with a second class particle is established and 
its relation to the KPZ pl| and noisy Burgers [E2[] equation are discussed. 



2 Model definition 

We consider two types of orders: limit orders that are wishes to buy or sell 
a given quantity of stock at a given price, and market orders that are orders 
to immediately sell or buy an asset at the best instantaneous price. Limit 
orders are stored in an order book until they are either canceled,[] or fulfilled, 
provided that the current market price has moved towards their prices. The 
model is constructed as a one dimensional lattice model, in a similar spirit 
as in |TJ|, |23]. Let the lattice of length L represent the price axis, with the 
lowest price on the left end at site 1 and the highest one on the right at 
site L. Limit orders of the two different kinds, i.e., asks (A) and bids (B) 
are placed on the lattice according to the price that the order is based on. 



^or instance because they had a predefined maximal lifetime. 
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As bids name lower prices than asks, they will be found on the left side of 
the lattice. The current market price (x) separates the two regions. Sites 
representing prices at which currently no order is placed are indicated as 0. 
In the model we made the following simplifying assumptions: 

• Only one kind of asset is traded and its price dynamics is not directly 
influenced by outside sources but just by the state of its limit order 
book. 

• Each site can only carry one order (exclusion model). 

• Limit orders of either kind come in a unit size. 

• Only a finite price interval of width L is considered. 

The last assumption is justified as trade only takes place in a narrow interval 
around the market price. In our model L can be chosen arbitrarily. For a 
discussion on differences between models and real limit order markets see e.g. 
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Figure 1: Example of a configuration and possible moves with their assigned 
rates. 



The dynamics of the lattice is as follows (see Figure 1): 

• At site 1 bids enter the system at rate a: — > B. 

• At site L asks enter the system at rate 7. — > A. 

• Asks and bids can diffuse one site towards the market price at rate 
p provided no other order is already placed at the target site: BO — > 
OB, OA -> AO. 

• Bids can be placed at unoccupied sites left of the market price at rate 
rj. The same holds for asks being placed right of the market price: 
-> A, -> B. 
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• Bids and asks can be evaporated at rate rj: A — > 0, B — > 0. This 
reflects both orders being canceled as well as being timed out. We 
have made the simplifying assumption of a constant removal rate at 
each site instead of considering the lifetimes of individual orders. 

• An ask can be fulfilled at rate p by an incoming market order, provided 
it is adjacent to the current market price. Thus the order is removed: 
xA — > Ox. 

• A bid can be fulfilled at the same condition and rate, leading to a 
decrease of the market price and removal of the order: Bx — > xO. 

The role of order injection and diffusion is to ensure a fluctuating order 
density on both sides of the price. On the other hand, the dynamics of 
the special particle which represents the price, is such that the sign of price 
increments is constant as long as the bid-ask spread is not minimal, i.e., as 
long as the price is not surrounded by two orders. This is a crude but efficient 
way of implementing trends in limit order market models. Notice that this 
dynamics implies that the price is always between the best bid and best ask 
orders, which is true 95% of the time in ISLAND ECN (www.island.com). 



Even if it is likely that orders do not diffuse [16], we use this ingredient as a 



way of obtaining exact results for the Hurst exponent. 



3 Simulation results 

Throughout our simulations we chose initial configurations where each site of 
the lattice is randomly occupied by an order with probability 1/2. Further- 
more we chose a = 7 = 1/4, p — 1/2. The lattice size L was chosen such that 
the market price M could not fluctuate out of the represented price interval 
during the simulation time. The choice of rates guarantees that the price has 
no drift but on average remains on its initial position, i.e., L/2. Averaging 
over initial conditions is implied in all our simulations. 

In the case rj = we remain with a model where price dynamics is solely 
caused by diffusion of limit orders. We are mainly interested in the Hurst 
exponent H defined by the relationship ((Sx) 2 ) 1 ^ 2 oc At H . 

In Fig. 2 we show the fluctuation of the price position x relative to the 
initial price versus time, in a double-logarithmic plot. As can be seen the 
Hurst exponent of the model is H = 2/3 for all times without a crossover. 
This behavior is in contrast to the corresponding result from the basic Bak, 



Paczuski, Shubik (BPS) model |14], 23|. In the BPS models offers to buy 



and sell diffuse on a lattice representing prices just as in our model. The 
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Figure 2: Fluctuations of the price sig nal y/(x 2 (t)) oc t H at 7] = and fit 
with H = 2/3. 



difference is that upon meeting offers to sell and buy mutually annihilate 
in the BPS model, thus carrying the model to the universality class of the 
A + B — > reaction diffusion model. For that model it is known analytically 
that H = 1/4 at long times plus logarithmic corrections at shorter times 
In our model no mutual annihilation (of ask and bid), takes place, but 
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just one type of order vanishes (ask or bid, fulfilled together with a market 
order), thus causing a price change. This carries our model into the realm of 
the KPZ universality class as we will illustrate in the next section and yields 
H = 2/3. 




Figure 3: Correlations of the increments at rj = and fit by the eye with a 
function decaying algebraically as t~ l l 2 . 



Price increments Sx(t) = x(t' + t) — x(t') show algebraically decaying 
correlations (Figure 3), with (5x(t f ) 5x(t' +t)) oc t -1 / 2 whereas these corre- 
lations should be essentially zero; this is due to the absence of evaporation 
(see below); the correlation of absolute increments has algebraically decreas- 
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ing autocorrelation with an exponent of approximately 1. These long ranged 
correlations cause the over- diffusive behavior. 
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Figure 4: Distribution of the increments (logarithmic scale) for t = 50 at 
7] = and Gaussian fit. The inset shows the variance of the distribution as 
a function of time and a fit proportional to t 4//3 . 

Due to the price process itself (|&r(l)| G 0,1), the histogram of 5x(t) 
is almost Gaussian in shape, the tails appear even less pronounced than a 
Gaussian pdf (Figure 4). The variance of the distribution of increments Sx(t) 
increases as a a (see inset of Figure 4). The dynamical exponent of the 
price process extracted from this property is z = 3/2. Clearly the stochastic 
process causing the price movements is not Gaussian, not even a rescaled 
Gaussian. 
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Figure 5: Fluctuations of the price signal y/ (x 2 (t)} oc t H at rj = 1/512 and 
fit with H = 2/3 in the early time regime and H = 1/2 at late times. 

In what follows we consider the case rj ^ 0. Clearly this is more realistic 
than rj = as it is possible to place orders of either kind at any unoccupied 
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site on the price axis without having to perform diffusion steps all the way 
from the boundaries. Also the withdrawal of orders due to cancellation and 
timeout becomes thus possible. As seen in Figure 5 the fluctuations of the 
price signal show a crossover from super- diffusive behavior at short times, 
characterized by the Hurst exponent H = 2/3 to diffusive behavior at later 
times, implying H = 1/2. 
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Figure 6: Decay of the correlations of the increments for 77 = (upper curve) 
and 77 = 1/8 (lower curve). 

The local random events controlled by the parameter 77 destroy the long 
ranged time correlations of the price increments. This can be seen in Figure 
6, showing the correlation of increments (5x(t) 8x(0)) as a function of t for 
77 = and 77 = 1/8. At t — 200 the correlation function for the 77 = case is 
still different from zero, while a decay to zero for the other case occurred long 
since. Note that the autocorrelation of price increments should be negative 
for short time whereas it is positive in our model; this is due to the fact that 
we do not allow the coexistence of the two types of markets orders (see [|30|1). 

Adjusting the rate 77 serves as a parameter to control the crossover point. 
Figure 7 shows the price fluctuations for values of 77 between 1/8 and 1/512. 
The larger the rate 77 for local events, the shorter is the time span for over- 
diffusive fluctuations. In fact, in our simulations over-diffusive behavior over 
a longer time interval appears only to be possible at seemingly meaningless 
low rates 77, namely 77 ~ 1/L, compared toa = 7 = l/4 and p = 1/2. From 
the study of empirical data of the Island ECN limit order market conducted 
in [[U| it is known that about 80 per cent of limit orders in the respective 
market vanished due to timeouts. Only 20 per cent of the offers were (at least 
partially) fulfilled. We measured these quantities as a function of 77 in our 
simulations, where fulfillment of an order means either the process Bx — > xO 
or xA — > Ox and timeouts are reflected by the rate 77. It turns out that for a 
lattice of L = 1000 at 77 = 1/512 about 8 per cent of orders were fulfilled and 
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Figure 7: Fluctuations y/(x 2 (t)) for 77 = 1/8, rj = 1/32, 77 = 1/128 and 
i] = 1/512 (from below). The parameter 77 controls the range over which 
H = 2/3 (dashed line). 

at 77 = 1/1024 about 16 per cent. Thus the choice of small 77 matching the 
observed fulfillment rate is realistic and could in fact be used as a means to 
gauge the simulation time by comparing the known empirical crossover time 
and the simulation crossover. 

In the spirit of dynamical scaling it is tempting to assume that the price 
fluctuations with sufficiently low 77 can be described in terms of a scaling 
function F with 

< x 2 (t) >= rf^Fi^t) 

Since 77 is a rate with inverse dimension of time one expects /i 2 = 1 for 
covariance of < x 2 (t) > under rescaling of time. For small times, i.e., small 
arguments of the scaling function, this ansatz should reproduce the behavior 
< x 2 {t) >oc t 4 / 3 which implies F(y) oc t/ 4 / 3 for y — » 0. Independence of 77 
thus yields the scaling relation 

4/3/i 2 = (1) 

Hence one expects ^ = 4/3. For large times crossover to diffusion implies 
F(%f) oc y for y — ► 00. For large 77 (of order 1 and larger) we obtain ordinary 
random walk behavior even at early times and the scaling relations are not 
expected to be valid. 

These arguments are well born out by Monte Carlo simulations. The best 
fit for the data could be achieved for the choice pix = 1.31 and /j 2 = 1-04. 
These exponents are used in the plot (Fig. 8). The scaling property suggests 
that the relevant time scale of the model is r = I/77, which is the average 
time between successive placements or evaporations of an order at a given 
site. 
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Figure 8: Data collapse for the data for 77 = 1/32, to r\ = 1/1024 using the 
scaling function given in the text and lines to guide the eye corresponding to 
H = 2/3 and H = 1/2. 
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Figure 9: Distribution of the increments for t=50 at rj = 1/8 and Gaussian 
fit. The inset shows the variance of the distribution as a function of time 
and the linear asymptote. 



The increments of the price signal 5x(t) for 77 7^ (Fig. 9) differ from rj = 
in two important respects. Firstly, the tails of the distribution are closer to 
a Gaussian. A second and more pronounced difference is the behavior of the 
variance a of the distribution as a function of time, which shows a crossover 
from a oc t 4//3 to a oc t just as the price signal itself. This means that the 
price performs an ordinary random walk at long times. A snapshot of the 
price movement is shown in Figure 10. 
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Figure 10: Snapshot of price movement at 77 = 1/8. 



4 Connection to the TASEP 



The virtue of our model is the equivalence at rj — to the totally asymmetric 
exclusion process (TASEP) for which a wealth of exact results exists P0fl . In 
the TASEP excluding particles enter a lattice at rate a from the left and 
hop with rate p to the right, provided the target site is empty. At the right 
end they can leave the system with rate (3. In connection with the TASEP a 
'second class' particle has been defined to have the following dynamics: 
A first class particle meeting upon a second class particle to its right will 
exchange places. A second class particle with a vacant right neighboring site 
hops to that site. The second class particles motion is designed such that 
it does not interfere with the motion of the first class particles. In fact, the 
motion of a single second class particle in the system is on average that of a 
density fluctuation in the system. 

Upon coarse graining the dynamics of the TASEP can be described by 
the noisy Burgers equation, which is closely related to the KPZ equation 
known to have a universal dynamical exponent 2 = 3/2 ||19|| . This implies a 
Hurst exponent H = 2/3 as discussed above. For the noisy Burgers equation 
the over- diffusive spreading of a density fluctuation (i.e., the spreading of 
the second class particle, representing the price signal in our model) with 
H = 2/3 has been shown analytically [2^] in the case of statistically averaging 
over initial positions as well as realizations, which is always implied in our 
simulations. 

The mapping between our model and the TASEP is as follows: The 
market price x represents the second class particle. Left of its position, bids 
are first class particles in the sense of the TASEP and vacancies remain what 
they are. To the right of the market price vacancies take the role of first class 
particles in the TASEP sense and asks that of vacancies. The price dynamics 
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in our model is precisely that of a second class particle or density fluctuation 
in the TASEP. 

The TASEP may be seen as a discretized version of the noisy Burgers 
equation. It is an exactly solvable model for which z = 3/2 has been ob- 



tained through the Bethe ansatz [27]. More recently, also the distribution 



of the second class particle for averaged random initial conditions has been 
calculated exactly through a correspondence with statistical properties of 
random matrices P5[ . This confirms the results of our simulations for a fi- 
nite lattice with open boundaries, but system large enough to be equivalent 
to an infinite system. We have also performed simulations for a fixed ran- 
dom initial state. We found that the super-diffusive spreading of the second 
class particle prevails, but the amplitude < x 2 > / 1 4 / 3 depends on the initial 
condition. This is in accordance with expectations |29[ . 



5 Conclusions 

In this work we have presented a model exhibiting the empirically observed 
crossover of the Hurst exponent from H > 1/2 to if = 1/2. By a mapping to 
the totally asymmetric exclusion process we obtain the exact value H = 2/3 
which is close to what is often observed in real markets. The existence of 
an exact analytic solution puts our model in contrast to the model by Bak 
I3| exhibiting over-diffusive spreading by volatility feedback into the system 



and a copying strategy of the traders, but for which no analytical solution is 
known. The over- diffusive behavior results from time correlations build up in 
the biased internal motion of asks and bids respectively which together with 
market orders drive the price process. We identify the average time between 
evaporation events of orders (due to time-out or cancellation respectively) 
at a given site as the relevant time scale of the model before crossover to 
diffusive Gaussian behavior. 
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